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Abstract. NGC 2359 is a Wolf-Rayet (W-R) nebula partially bound by a rather dense and warm molecular cloud. We present 
the results derived from CO and 13 CO fully sampled maps of the molecular material with angular resolutions up to 12". We have 
detected three different velocity components, and determined their spatial distribution and physical properties. The kinematics, 
morphology, mass and density are clearly stratified with respect to the W-R star. These features allow us to learn about the recent 
evolutionary history of HD 56925, because the multiple layers could be associated to several energetic events which have acted 
upon the surrounding circumstellar medium. Hence, a careful study of the different shockfronts contain clues in determining 
the present and past interaction of this evolved massive star with its surroundings. From the analysis of the mass-loss history in 
massive stars like HD 56925, we suggest that the multiple layers of shocked molecular gas are likely to be produced during the 
earlier LBV phase and/or the actual W-R stage of HD 56925. 
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1. Introduction 

NGC 2359 is one of the archetypal wind-driven Wolf-Rayet 
(W-R) nebula (Chu et al. 1983). It is an optical neb- 
ula excited by the WN4 star HD 56925 ((RA, Dec.) 2 ooo = 
(OVMS^.O.-Ona'Ol'.'S)), which is catalogued as WR7 
by van der Hucht (2001). NGC 2359 has an almost spherical 
morphology, and contains some N-enriched filaments inside 
(Esteban et al. 1990). The southern part of the nebula -which 
contains most of the associated molecular gas- has been ob- 
served in the recent years. 

St-Louis et al. (1998) have detected the H 2 1-0 S(l) emis- 
sion in coincidence with Ha, but they could not disentangle 
the dominant excitation mechanism of the molecular hydrogen. 
Carbon monoxide was firstly detected by Schneps et al. (1981), 
and subsequently mapped in the J = 1 — » and 2 — > 1 lines by 
Rizzo et al. (2001a, Paper I) and by Cappa et al. (2001). While 
Cappa et al. focussed in determining the total neutral mass and 
in the energetics, Paper I has concentrated on the interaction 
processes between the nebula and the surrounding molecular 
gas. 

In Paper I, it was shown that the ambient molecular cloud, 
observed at a velocity in the Local Standard of Rest (Vlsr) of 
67 km s _I , was accelerated to Vlsr = 53-54 km s by the stel- 
lar wind from HD 56925. It was also established the presence of 
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a rather hot region just coincident with the H 2 vibrationally ex- 
cited emission surrounding the southern optical nebula, as well 
as a high-opacity zone adjacent to the cold quiescent cloud. 
Observations of 13 CO have revealed another striking feature: 
in the region of higher intensity, the 13 CO J = 1 — > line peaks 
at 54-55 km s" 1 , 1-2 km s" 1 higher than the CO J = 1 -> 
and the J = 2 — > 1 lines. It is worth mentioning that, at the 
same position, the metastable (1,1) and (2,2) lines of ammonia 
(NH3) were detected for the first time in a W-R nebula (Rizzo 
et al. 2001b). From the NH3 data, a rotational temperature of 
30 K and a relatively high abundance of 10~ 8 were determined. 

These observational findings indicate the nature of the in- 
teraction with the ambient material during different evolutive 
stages from the hot star. The molecular cloud has been exposed 
to the UV radiation from the hot star and to the shocks pro- 
duced by the stellar winds from the W-R star and the mass ejec- 
tions from previous stages. According to the hydrodynamical 
models of Garcfa-Segura et al. (1996a, 1996b), several shocks 
are expected at a few pc from a massive evolved star, as a 
consequence of the time-dependent stellar winds of the post- 
main sequence evolution, mainly the red supergiant (RSG) and 
the luminous blue variable (LBV) phases. Hence, the molecu- 
lar cloud surrounding NGC 2359 provides us with an excellent 
laboratory to study the interaction between the massive stars 
and their environment. W-R stars will disturb the circumstel- 
lar medium (CSM), generating an ionized nebula surrounded 
by a shocked region. In particular, the NH3 would be a good 
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molecule to disentangle the effect of the radiative field and 
the shock in the environment of W-R stars. The relatively high 
abundance of NH3 strongly suggests a shock chemistry in the 
outer part of the cloud, since it would be released from the dust 
grain mantles by non-dissociative shocks (Flower & Pineau des 
Forets 1994; Rizzo et al. 2001b). 

A detailed knowledge of the spatial and kinematical distri- 
butions of the CO isotopes would allow us to know with more 
precision the physical conditions and the kinematical features 
in this scenario, together with an approach to the knowledge of 
the dominant mechanisms of excitation (UV fields or shocks). 
In order to satisfy these goals, we combine in this paper high 
angular resolution CO J = 1 -» 0, CO J = 2 -> 1 and 13 CO 
J = 1 -> maps with high sensitivity spectra of CO, 13 CO 
and C ls O taken at selected positions. This observational study 
reveals the presence of hot dense regions, warm regions and, 
more conclusively, the stratified distribution of different kine- 
matical components, associated to several shockfronts which 
have been acting in the recent past of HD 56925. 



2. Observations 

We have used the IRAM 30-m radiotelescope at Pico Veleta 
(Spain), during two observing runs in June and August, 2001. 
We have performed two different observing modes, namely on- 
the-fly maps and single position observations. 

The on-the-fly maps were obtained for the CO and 13 CO 
J = 1 — > 0, and the CO J = 2 — » 1 lines simultaneously. We 
have mapped a 240" x 120" area toward the southern part 
of NGC 2359. The scanning directions were orthogonal in the 
equatorial system, e.g., at constant right ascension and constant 
declination. Row spacing between scans were 6", and there- 
fore all spectral lines maps are fully sampled. Individual maps 
were later combined using the GILDAS package. One receiver 
was tuned for each J = 1 — > line, while the CO J = 2 — > 1 
line was observed simultaneously by two different receivers. A 
256-channels filter bank with 100 kHz of channel resolution 
and a 2048-channel autocorrelator were used as backends, pro- 
viding velocity resolutions between 0. 1 and 0.3 km s _1 . Typical 
rms was 0.3 K for all the observed lines. All the temperatures 
throughout this paper are referred to an atmosphere corrected 
scale, e.g., T^. 

The single position observations were performed along a 
strip of particular interest in position-switching mode. We ob- 
served the J = 1 -> and J = 2 -> 1 lines of CO and 13 CO 
simultaneously, with integration time between 15 and 40 min- 
utes. As a byproduct, the C ls O J = 1 -> line was also ob- 
served with the same receiver than the 13 CO J = 1 — > line. We 
used the same backend as for the on-the-fly maps, and hence we 
had identical velocity resolutions. The achieved rms was better 
than 50 mK in all cases. 

Calibration was performed during the observations and 
later checked up on standard sources (Mauersberger et al. 
1989). The HPBW was 22" for the CO J = 1 -> line, 11" 
for the CO J = 2 -» 1 line, 24" for the 13 CO J = 1 -> line and 
12" for the 13 CO J = 2 -> 1 line. 



3. Results 

3.1. On-the-fly maps 

The southern part of the nebula was mapped in the CO J = 
1 -> 0, J = 2 -> 1 and 13 CO J = 1 -> lines simultane- 
ously. In Paper I, it has been shown that this region contains 
the most intense and interesting CO features, between 52 and 
67 km s _1 . The maps in Fig. 1 show the integrated emission of 
these lines from 50 to 56 km s _1 . Figure la shows the whole CO 
J = 2 — > 1 emission from Paper I, superimpossed to an optical 
image of the nebula. The same figure also indicates the region 
mapped in this work. Figures lb-Id show the results of these 
data in the three spectral lines already mentioned. We can see 
that the molecular cloud perfectly matchs the southern border 
of the optical nebula. The CO and 13 CO emission are confined 
to layers which are less than 70" wide. 

The overall emission maps depicted in Fig. 1 do not 
show significant differences among the observed spectral lines. 
However, a number of new features appear when the emission 
is analysed in smaller velocity bins. Figure 2 shows an overlay 
of the CO J = 1 — > and J = 2 — > 1 emissions in two different 
channel at velocities of 55.2 km s _1 (grey contours) and 50.1 
km s _1 (black contours). It is clearly noted a separation of more 
than 20" between the emission at both velocities. Actually, the 
thinner layer located toward the inner part of the nebula, at a 
velocity of ~ 50 km s _1 , is adjacent to the more extended emis- 
sion. The width of this layer is of about 15 - 20", as we can see 
in the CO J = 2 — > 1 map. 

To better understand the velocity structure of the CO emis- 
sion we have analyzed the data along several position-velocity 
strips. The most illustrative one is shown in Fig. 3. This strip 
crosses both the large and the thin molecular layers. It is clearly 
shown a sudden line broadening just at the position of the thin 
layer, around 50 km s _1 . This behaviour clearly demonstrates 
that the feature noted in Fig. 2 is not merely due to a velocity 
gradient of a single component, but to the appearing of a second 
component only present at the interface between the molecular 
cloud and the optical nebula. 

3.2. Individual positions 

In order to have a complete picture of this striking feature, we 
performed high-sensitivity, individual observations at the po- 
sitions shown as squares in Fig. 3a. These positions, labelled 
with letters from A to E, are spaced by ~ 20". Observed spec- 
tral lines at these locations are the CO J = 1 — > 0, CO J = 2 — > 1, 
13 CO J=l -> 13 CO J = 2 -> 1, andC ls O J=l -> lines. 
The 13 CO J = 2 -> 1 and C ls O J = 1 -> lines were detected 
for the first time in a W-R environment. 

The resulting spectra are shown in Fig. 4. For all the ob- 
served lines, the peak velocity remains roughly constant from 
positions A to C, and the linewidth increases. Both CO lines at 
position C are stronger than at position B, while in both 13 CO 
lines we note exactly the contrary. Toward position D all the 
line intensities are roughly reduced to one half and a strong as- 
symetry and line broadening appear. There is significant emis- 
sion at velocities as low as 47 km s _1 . The broadening shown in 



Rizzo et al.: Molecular shocks in NGC 2359 



3 



the position-velocity strip (Fig. 3b) is coincident with this posi- 
tion. C ls O is only detected in positions B and C, corresponding 
to the highest emission at the other isotopes. 

Position E, although weak, shows a large CO J = 2 — > 1 
broadening to velocities as low as 42 km s _1 . This implies a 
total velocity interval emission of at least 14 km s _1 . It is worth 
mentioning that this very broad emission at V L sr < 50 km s _1 
is not detected in the CO J = 1 — > line. As mentioned above, 
we see in Fig. 3b a notable line broadening between positions 
C and D, where the lower levels are seen at 47 and 57 km s _1 . 
Surprisingly, the high sensitivity individual positions reveals 
a third velocity component at position E, which could not be 
detected using the on-the-fly maps because of the sensitivity. 

The main observational findings are summarized in Table 
1 . For each position, this table shows the parameters of the five 
observed transitions, obtained from a gaussian fit to a single 
component: antenna temperature of the maximum (T p k), LSR 
peak velocity (Vlsr) and full width at half maximum (AV). 
The single component assumption is obviously not valid in 
some cases due to the strong assymetries of the line profiles, 
and hence is only provided for comparison. 

4. Physical parameters 

4.1. Velocity components 

The observational features shown in Sect. 3 (in particular 
Figs. 3 and 4, and Table 1) lead us to define three different 
velocity components. The first component corresponds to the 
most extended molecular gas and is the broad component, al- 
ready studied in Paper I, at velocities from 52 to 57 km s _1 . 
In the following, we shall refer to this component as the "low 
velocity component" (Ivc). The higher spatial resolution of our 
maps allow us to disclose the presence of a second, less ex- 
tended velocity component, seen as a thin layer (~ 20" width) 
in Fig. 2, in the velocity range 48-52 km s _1 . We shall refer to 
this component as the "intermediate velocity component" (ivc). 
Finally, the high sensitivity of the individual spectra revealed 
the third component, hereafter referred to as the "high veloc- 
ity component" (hvc), as a very compact one, clearly detected 
only in the CO J = 2 — > 1 line toward position E (Fig. 4), at 
velocities lower than 48 km s - ' . 

The extension in the sky corresponding to each component 
may be roughly inferred from Fig. 3(a): the Ivc is sketched in 
the greyscale, the ivc is plotted in contours, and the hvc is only 
detected in position E. The Ivc corresponds to most of the gas 
detected in Paper I. It is the most extended component, and 
bounds the whole optical nebula by its southern and eastern 
borders. The ivc is at the interface between NGC 2359 and the 
most intense emission of Ivc. It seems to form a narrow layer 
between the Ivc and the optical nebula. Finally, the position 
where the hvc was detected is projected onto the ionized gas 
and presumably has a compact morphology. 

The stratified location of the three components with respect 
to the optical nebula NGC 2359 and the quiescent molecular 
gas at 67 km s _1 , might be explained by succesive events as- 
sociated to the evolution of a massive star. Therefore, different 
ages are expected for each component. The Ivc is the most ex- 



tended, uniform and massive one, its velocity is the closest to 
the quiescent gas and hence is thought as the oldest one. The 
ivc, however, is less extended than the Ivc and closer to the 
W-R star. These kinematical and morphological features might 
be explained by an event faster and younger than that which 
formed the Ivc. So far, the ivc may be the result of the encounter 
of two winds, in a manner suggested by Garcfa-Segura et al. 
(1996a, 1996b). Finally, the hvc may be recently shocked and 
would be the youngest component. 

4.2. Overall features 

Our on-the-fly CO and 13 CO maps trace the distribution of the 
molecular gas both into and around the southern optical nebula, 
including the interesting molecular interface at its border. The 
different velocity components have probably been produced by 
violent energetic mechanisms usually linked to the massive star 
evolution. During the evolution of HD 56925, the stellar wind, 
eventual envelope ejections and the strong UV radiation fields 
could have produced shocked, ionized and photodissociated re- 
gions. The relative variations of different spectral lines should 
tell us about the physical conditions of the molecular gas. 

We have analyzed the line ratios among the three observed 
spectral lines which covers the southern and eastern borders 
of NGC 2359. In the following, we shall refer to the line ratio 
of 12 CO J = 1 -> to 13 CO J = 1 -> as R isot , and to the 
line ratio of 12 CO J = 2 -> 1 to 12 CO J = 1 -> as R exc . 
Both R; sot and R exc were computed for different velocity ranges 
and for positions with line intensity greater than 3<x. R exc was 
computed after convolving the CO J = 2 — » 1 map to the J = 
1 — > angular resolution. While R exc is a good tracer of the 
excitation temperature, R; sot might trace the combined effect of 
isotopic variations and changes in the CO opacity. 

Figures 5a and 5b show the spatial distribution of R exc for 
the velocity ranges (52, 55) and (48, 51) km s _1 , roughly equiv- 
alent to the Ivc and the ivc, respectively. For the Ivc (Fig. 5a), 
R exc does not significantly changes in the brightest zone, and 
it is probably greater toward the south. A mean value for R exc 
of 0.9 and 1.2 is found to the eastern and southern part of the 
map, respectively. We have obtained similar results for the ivc 
(Fig. 5b), where a mean value of 1 is found. However, the 
13 CO emission is really striking when compared to the CO J 
= 1 — > emission. Figure 6 sketches the distribution of R; sot in 
the velocity interval (52, 55) km s _1 (roughly the Ivc). It clearly 
changes sistematically from southeast to northwest, and varies 
in almost one order of magnitude. The sensitivity of our on- 
the-fly maps were not enough to detect the 13 CO emission in 
the velocity range corresponding to the ivc, and R; sot could not 
be computed. 

4.3. Excitation conditions 

From the overall information shown in Figs. 5 and 6, we see 
a rather uniform excitation conditions over the whole region. 
So far, in most of the observed area, the excitation temper- 
ature does not varies significantly. However, large variations 
in the CO opacity might be inferred. If we assume a constant 
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CO/CO isotopic abundance, the CO lines will be optically 
thick in the eastern part of the cloud. This feature explains the 
velocity difference between both isotopes, as noted in Paper I. 

We have estimated the physical conditions for the three ve- 
locity components, by modelling the emission of the observed 
lines for positions A to E, using the large velocity gradient 
method (LVG). The LVG method implements the Sobolev's 
formulation developed by Castor (1970) and later extended by 
Jeffery (1995), among others. We have computed the parame- 
ters assuming kinetic temperatures (T K ) from 10 to 100K. For 
all components but the live, a minimum T K of 30 K is needed 
to match the observed lines. This lower limit for Tk is in good 
agreement with the results of Rizzo et al. (2001b), based on 
the metastable NH 3 (1,1) and (2,2) lines detected at position B. 
Furthermore, the results are not critically different for the tem- 
perature range 30-100 K. Both CO lines result to be optically 
thick in most cases, whereas both 13 CO lines were optically 
thin in all cases. In consequence, we used the 13 CO 2-1/1-0 
line ratio (on a main-beam temperature scale) to estimate the 
H2 density, n(H2). For a given T K , the column densities of 13 CO 
and C ls O were determined from the derived n(H2) and the J = 
1 — > line intensity of each isotope. The CO column density 
was estimated by assuming a constant isotopic ratio of 60 be- 
tween CO and 13 CO (Wilson & Matteucci 1992) 

The Table 2 shows the results of applying the LVG method 
to each component at each position, by assuming a T K of 30 K. 
By adopting this value for Tk we are in the most conservative 
case, because T K may be higher in position E and probably 
in position D. Anyway, the column densities do not vary in 
more than 10 percent (significantly smaller than other sources 
of uncertainty) in this temperature range. The first and sec- 
ond columns indicate the position and velocity component. The 
third column provides the obtained n(H2). From fourth to sixth 
columns, column density of CO, 13 CO and C ls O are provided, 
respectively. 

The results derived for hvc in position E need to be ex- 
plained in detail. This velocity component was only detected in 
the CO J = 2 — > 1 line, and therefore has a minimum R exc of 6. 
Such a high value could not be obtained having a Tk of 30 K, 
even by considering that the emiting region is unresolved by 
the 12" beam of the CO 2 — > 1 line. Therefore, it is neccesary 
larger values of T K and a higher density to match the observed 
ratio. We estimate that T K must be larger than 80 K for this 
component, as well as the density should be at least 10 5 cm 4 . 
The undetection of this component in the CO J = 1 — > line 
strenghtens the idea of a small emiting volume, and a higher ex- 
citation temperature is also expected in this position (towards 
the inner part of the ionized nebula), as discussed in Paper I. 
Finally, it should be noted that the CO column density for this 
component has been estimated from the CO J = 2 — > 1 line. 

5. Origin and evolution of the CSM around WR7 

5.1. Dynamics of the molecular gas around NGC 2359 

The morphology and kinematics of each component, as well 
as their location with respect to NGC 2359 and the quiescent 
molecular gas have been discussed in previous sections. Some 



of the physical conditions have been summarized in Table 2. 
In this section, we will put them together in order to draw a 
picture of the region and its possible evolution. 

In Fig. 7, we plot the variation of the gas density and 13 CO 
column density as a function of the projected distance to the 
W-R star HD 56925, assumed at a distance of 5 kpc (Goudis et 
al. 1994). The plots correspond just to the Ivc and the ivc and 
not to the hvc, because it is only detected at a single position. 
The stratification of these components are not only evident in 
the morphology and the kinematics. This progressive variation 
is also present in the maxium density (roughly 10 3 crrT 3 at Ivc, 
10 4 cirT 3 at ivc and at least 10 5 cirT 3 at hvc) and in the total 
mass involved (regarding the column densities and the spatial 
extension of each component). If we take into account the ob- 
served low values of R; sot , it is also expected a trend in the opac- 
ity, decreasing from Ivc to hvc. The temperature is expected to 
increase toward the W-R star, and therefore in the sense Ivc- 
ivc-hve. For the hvc, both the density and the temperature must 
be higher than in the other two components. Although a care- 
ful chemical study of the region is needed, we can associate 
each component to different evolutive episodes previous to the 
actual W-R status of HD 56925. 

In Paper I, we have discovered a quiescent molecular 
cloud, that was interpreted as the original gas, which was later 
shocked by the expansion of a wind previous to W-R. The pres- 
ence of molecular gas close to a W-R star is indeed puzzling, 
because it is located at a projected distance of only ~ 8 pc, and 
has been exposed to the very strong radiation and stellar winds 
during all the evolutive phases of HD 56925. Therefore, the 
three velocity components reported here are the result of such 
interaction and represent the disturbed molecular gas linked to 
some of those phases. In order to associate these features with 
the W-R phase or a previous one, we have made crude esti- 
mates of some global dynamical parameters for each compo- 
nent, which are provided in Table 3. The whole nebula and its 
surroundings was completely mapped in Paper I, and we con- 
sidered that the mass obtained there includes both the Ivc and 
the ivc. In this work, we roughly estimate that 80% of that mass 
are from the Ivc, while the other 20% are from the ivc. The mass 
of the hvc was estimated from the CO J = 2 — > 1 spectra and 
assuming an angular extension of one beam. The computation 
of masses and projected distance were done by assuming a dis- 
tance of 5 kpc (Goudis et al. 1994). The expanding velocities 
were derived by assuming a rest velocity of 67 km s _1 (Paper I). 
The dynamical time was computed according to Dyson (1989). 
Finally, Table 3 also includes estimates of kinetic energy and 
momentum of each component. 

The mass is significantly different for each component, de- 
creasing in the sense Ivc-ivc-hvc, while the other parameters 
are of the same order of magnitude. All the components have 
dynamical ages of at most 10 s yr. This time represents an up- 
per limit of the age of the shock which produced the actual 
kinematics. It is one order of magnitude lower than the main- 
sequence Hi shell (2.3 10 6 , Paper I), and quite larger than the 
dynamical age of NGC 2359, 1.3 1 5 , determined by Treffers 
& Chu (1982). However, it should be noted that this dynami- 
cal age for the optical nebula is indeed an upper limit, which 
is reflected by the 1.3-2.2 10 4 yr determined from Hipparcos 
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proper motion for HD 56925 between the center of the nebula 
and its actual position. Due to the high symmetry of NGC 2359, 
we think that a few 10 4 yr represents a reasonable approach of 
the actual NGC 2359's age. In consequence, the observed CO 
features have been produced at an epoch between the end of the 
main sequence (MS) and the beginning of the W-R phase. 

According to the current models, a W-R is the descen- 
dant of an O-star of at least 30 M . Depending on its initial 
mass, the O-star may evolve to W-R through different evo- 
lutive paths, including a red supergiant (RSG) or a luminous 
blue variable (LBV) phase. Other intermediate stages may ac- 
count, such as Of or different sub-classes of W-R (WNE, WNL, 
WC). Although the MS phase is indeed a long-lived one, char- 
acterized by a strong stellar wind and an enormous rate of UV- 
radiation, we could not simply relate the features observed in 
the molecular gas around NGC 2359 to this phase, untill the 
other phases were also considered. Moreover, a MS bubble 
around HD 56925 was found in Hi (Paper I), and it has a lin- 
eal size of 35-40 pc, closely similar to many other cases in the 
Galaxy, and in good agreement with the theory. 

It is indeed possible that the interaction of the molecular 
cloud with the MS wind could not be strong enough and the 
stellar wind could ran through the molecular cloud. This is 
a reasonable argument owing to the relatively low density of 
the MS wind. We now consider the other evolutionary phases 
which may account for the observed features. The more likely 
mechanism is that the stellar wind has accelerated and heated 
the quiescent molecular cloud, and hence it is more relevant the 
instantaneous effect than the accumulated one along the whole 
phase. Therefore, the phase responsible of accelerating the gas 
should have a large momentum rate (P) and perhaps a large 
mechanical luminosity (L mec h). We also should pay attention 
to the density of the stellar wind, because it is an indicator of 
how large is, in comparison, the mass injected to the CSM at 
supersonic velocities. 

5.2. Models of stellar evolution 

Let us consider two different models of stellar evolution, which 
reproduce two standard evolutive paths: (1) the 35 M model 
of Garcfa-Segura et al. (1996b), with the path 0->RSG^W- 
R; and (2) the 60 M model of Langer et al. (1994), with the 
path 0->pre-LBV->LBV->W-R. We should note that some of 
these phases really includes several other, although no signifi- 
cant changes in the global stellar parameters are noted. So far, 
the W-R phase in both models can be divided in three sub- 
stages, such as a H-poor, a H-free and a WC stage. In model 
(2), the pre-LBV phase includes both the Of and the H-rich 
WN phase. The evolution of the CSM around both type of stars 
have been modelled by Garcfa-Segura et al. (1996a; 1996b). 

The Table 4 provides mean representative values from each 
phase in each model, such as the time spent, mass loss rate (M), 
L m ech, P> kinetic energy and total momentum (P) provided 
during the whole phase, and the wind density (p w ind) at 5 pc 
from the star. The 35 M model is split in two cases, having a 
slow (15 km s _1 ) and a fast (75 km s~') RSG wind. 



5.3. Acummulative effects 

In the three models, both the MS and the W-R stages provide 
the greatest fraction of the total energy injected to the envi- 
rons (10 50 ~ 51 erg). In the case of the MS stage, however, this is 
due to an accumulate effect, because the star spents some Myr 
in MS, one or two orders of magnitude larger than the other 
phases. In contrast, the W-R stage injects to the medium a sim- 
ilar amount of kinetic energy in some 10 5 yr, because its huge 
mass-loss rate. On the other hand, the total energy ouput asso- 
ciated to the RSG (fast wind) and LBV phases are two orders of 
magnitud (10 48 erg) lower than the MS or the W-R stages. The 
amount of kinetic energy associated to the RSG (slow wind) 
phase is even lower by about 4 orders of magnitude (10 46 erg). 

A similar situation occurs when considering the total mo- 
mentum provided in each phase, although the contrast is lower 
than in the energy. 

5.4. Instantaneous effects 

When we take into account the instantaneous effects (well de- 
scribed by L m ech, P and p w ind), the situation changes dramati- 
cally. The intermediate stages RSG and LBV are characterized 
by an increase of M by at least one order of magnitude with 
respect to their previous phases, and a sudden decrease of the 
terminal velocity. As a result, L me ch remains almost constant, 
while both P and p w ; n d increases significantly. 

During these short-lived phases, the circumstellar gas is 
filled with a large amount of material (~ 20 M ) from the stel- 
lar interior. This gas is relatively dense even at a few pc from 
the star and is moving at low (although still supersonic) veloc- 
ities. The RSG wind (slow wind) is the densest one, but the 
momentum rate is practically the same as in the MS case; so 
far, the effect onto the quiescent molecular cloud should not be 
different from that of the MS stage. The RSG wind (fast wind) 
is somewhat higher, but not too different from MS wind; a fac- 
tor of 2 or 3 in these simplified estimates should not be seen as 
a significant one. 

In contrast, the LBV wind is significantly different from 
the MS wind. The momentum rate is one order of magnitude 
greater than in any other evolutive stage, and the wind density 
is three orders of magnitude larger than the previous phases 
(MS and pre-LBV) and two orders of magnitude larger than 
the W-R stage. 

5.5. Interpretation 

If the progenitor of the W-R star HD 56925 has had a mass 
close to 35 M , the W-R wind is the unique source capable of 
disturbing the molecular cloud, because the RSG wind had not 
power enough to accelerate the original cloud. So far, the W- 
R wind had encountered an already stopped dense RSG wind 
in contact with the molecular cloud. Consequently, the W-R 
wind has been suddenly decelerated and produced the observed 
shocked molecular gas. 

If, on the other hand, the progenitor of HD 56925 has had 
a higher mass, it is very probable that the star had developed 
a LBV phase, which is characterized by a rather dense wind 
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and a high P. The LBV wind might produce what we are ob- 
serving in CO. The values provided in Table 4 are mean values 
for each phase as a whole. This simplification is good enough 
for all phases but LBV, because the parameters does not sig- 
nificantly change during the whole phase. However, the LBV 
phase consists of several episodes of enlarged stellar winds (or 
mass ejections) produced by hydrodynamical instabilities in 
the stellar structure. It is thought that the mass-loss rate during 
these episodes reaches values as high as several 10~ 3 M Q yr~', 
increasing the mean values of M, P, L mec h and p w i„d in Table 4 
by a factor 10-20 during some 10 3 yr. As in the 35 M case, 
the W-R wind may have found the RSG wind and produced 
the shock observed around the W-R nebula. The presence of 
several shockfronts (the three CO components and the ionized 
layers into the optical nebula) strenghten this idea. If only the 
W-R stage had produced the observed features, it should have 
similar episodes of increased mass loss. 

The whole optical nebula would be associated to the W- 
R phase, and its hot material (more or less uniform) would 
be located inside a multiple-layers wind with an earlier ori- 
gin. In this case, the bounding of the nebula can be atomic or 
molecular, depending on the previous presence of a molecular 
cloud close to the star. The morphology of the molecular gas 
in NGC 2359, as well as in NGC 6888 and the nebula around 
WR134 (Rizzo et al. 2003a) are strongly suggestive in this 
sense. 

In any case, it is clear the need of enlarge the sample, search 
for possible chemical differences which may confirm our hy- 
pothesis, observe and characterize the stellar winds of the well- 
known LBVs (P Cyg, rj Car, HR Car) and construct hydrody- 
namical models which take into account the different evolutive 
phases and the presence of a close molecular cloud. 

6. Future perspectives 

The history of a massive star is written by the impacts onto 
the surrounding gas, particularly in the molecular component. 
Therefore, we can increase our knowledge about massive stars 
evolution by observing the consequences in their environments. 
There are just a few observed cases of interaction between W-R 
stars and the surrounding molecular gas (Marston et al. 1999; 
Paper I; Cappa et al. 2001; Rizzo et al. 2003a). However, we 
think that NGC 2359 becames the best studied case, where a 
time evolution may be followed and where high-density gas 
has been detected. New observational efforts are neccesary in 
order to increase the sample and confirm some of the results 
obtained here and predicted by the theory. 

The increasing detection of molecular gas around W-R neb- 
ulae opens questions which are not completely fulfilled by cur- 
rent models. In the future, new hydrodynamical simulations 
should address this aspect of the interplay between evolved 
massive stars and their environs. It is expected to know where 
and when the molecules can be formed or can survive around 
these hot massive stars. 

On the other hand, the detection and chemical analysis of 
complex molecules in these ambients, at mm and sub-mm lines, 
would tell us about the chemical composition and the excitation 
conditions of the gas. In this sense, just a few detections have 



been reported (Marston 2001; Rizzo et al. 2001b, 2003b). This 
is a key subject, because we can learn not merely about the 
interplay of the massive star evolution with its surroundings, 
but also the ambient where a supernova takes place. So far, 
this field is also connected with the chemical evolution of the 
Galaxy and its energetic balance. 



7. Conclusions 

As shown in Paper I, the southeastern border of NGC 2359 con- 
tains the most clear sign of interaction of a W-R nebula with the 
molecular environment. The original molecular cloud was par- 
tially accelerated by at least 14 km s _1 (from 67 to 53 km s~') 
during one of the early evolutive phases of HD 56925. In this 
work, we have used the on-the-fly tecnique to obtain detailed 
maps of the CO and 13 CO emission of this disturbed molecular 
gas. The angular resolution and sensitivity achieved allowed us 
to distinguish three velocity components, clearly stratified with 
respect to the nebula and the quiescent gas. 

The first component (called the Ivc throughout the paper) is 
the most extended and the most massive one. The second com- 
ponent, the ivc, appears as a thin layer (less than 1 pc thickness) 
located at the interface between the Ivc and the optical nebula. 
Specially striking is the third component, the hvc, which ap- 
pears as a very compact component, only detected in the CO J 
= 2 — » 1 line towards the inner part of the optical nebula. The 
Ivc has LSR velocities between 52 and 57 km s _1 (the nearest 
to the quiescent gas), the ivc appears at 48-52 km s _1 , and the 
hvc emits at velocities between 42 and 48 km s _1 . So far, the 
stratification of these components are both morphological and 
kinematical features. 

By observing single positions with higher sensitivity along 
a strip -more or less radially with respect to WR 7-, some phys- 
ical parameters (volume density, column density and opacity) 
have been inferred. The density increases in the sense Ivc-ivc- 
hvc from 10 3 cirT 3 to at least 10 5 cm -3 , while the column den- 
sity, total mass and opacity should decrease in the same sense. 
We have proposed a scenario where the three components are 
produced during different evolutive episodes of the predeces- 
sors of HD 56925. Each component would be the result of an 
encounter (likely shocks) of the expanding stellar winds (or 
mass ejections) with the material ejected during the previous 
phase or even with the ambient gas. This idea qualitatively fits 
some theoretical results (Garcfa-Segura et al. 1996a, 1996b), 
although some specific aspects should be analyzed in detail. 
Particularly, the possible presence of molecular gas at a few pc 
from the massive star should be considered in future theoretical 
works. 
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Table 1. Gaussian fits to the observed CO isotopic spectral lines 



Pos. 


COJ = 


l->0 


COJ = 


2->l 


13 COJ = 


l->0 


13 COJ = 


2->l 


c 18 o J= 1^0 




T p k 


Vlsr 


AV 


Tpk 


Vlsr 


AV 


Tpk 


Vlsr 


AV 


Tpk 


Vlsr 


AV 


1 pk v LSR " v 




K 


km s 1 


K 


km s 1 


K 


km s 1 


K 


km s 1 


K km s -1 


A 


2.6 


53.4 


2.1 


2.6 


53.3 


2.3 


0.5 


53.8 


1.9 


0.4 


53.2 


2.2 


<0.04Kkms-' (fl) 


B 


6.3 


53.7 


2.8 


4.4 


53.6 


3.1 


1.5 


54.0 


2.1 


1.8 


53.7 


1.9 


0.08 53.6 2.6 


C 


7.4 


53.3 


4.2 


5.4 


53.0 


4.7 


1.0 


53.9 


2.7 


0.8 


53.2 


4.0 


0.08 53.7 2.8 


D 


3.6 


52.3 


5.3 


2.6 


51.8 


5.5 


0.3 


52.5 


5.8 


0.2 


52.0 


5.2 


<0.09Kkms-' (fl) 


E 


0.9 


53.2 


4.8 


0.6 


53.0 


14.0 


0.1 


53.7 


4.3 


0.1 


53.4 


1.5 


< 0.05 K km s" 1 (fl) 



<a) Upper limit of the integrated line intensity. 
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Table 2. LVG results for the CO isotopes in NGC 2359 



r OS. 


C-Omp. 


10 3 cm" 3 


10 16 cm" 2 


IN ^ v V ) ) 
10 14 cm" 2 


IN \J) 

10 14 cm" 2 




A 


Ivc 


3.0 


6.1 


10.2 


< 0.8 


> 12.9 


B 


Ivc 


4.9 


24.4 


40.7 


2.2 


18.2 


C 


Ivc 


3.0 


17.3 


28.8 


1.9 


14.8 


C 


ivc 


11.5 


3.8 


6.3 


< 1.5 


>4.3 


D 


Ivc 


3.9 


3.5 


5.9 


< 1.3 


>4.4 


D 


ivc 


3.0 


5.3 


8.9 


< 1.4 


> 6.5 


E 


Ivc 


1.4 


1.7 


2.8 


< 0.8 


> 3.5 


E 


ivc 


1.7 


0.2 


< 0.7 


< 1.0 




E 


live 


> 100 


0.2 


< 2.0 


< 1.6 
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Table 3. Global parameters of the molecular layers in 
NGC 2359 



Comp. 


Mass 




Rproj 


Tdyn 


E k 


P 




M 


km s -1 


pc 


yr 


erg 


M Q km s -1 


Ivc 


130 


12 


6.5 


3.010 s 


2.0 10 47 


1.710 3 


ivc 


30 


17 


5.8 


1.910 s 


1.210 47 


6.8 10 2 


hvc 


8 


23 


5.1 


1.210 5 


4.810 46 


2. no 2 
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Table 4. Evolutive phases of massive stars 



Phase Time M L mcch 


P 


E k 


P 


Pw'md 


yr M yr 1 erg 


M km s _1 yr -1 


erg 


ivIq Km s 


cm 


35 M G , slow wind. Model from Garcia- Segura et al. 1996b 
O-star 4.5 1 6 5.61(T 7 1.1 10 36 1.410~ 3 
RSG 2.310 5 8.110 5 5.810 33 1.210" 3 
W-R 1.910 s 2.410~ 5 3.410 37 5.010- 2 


1.610 s0 
4.210 46 
2.010 s0 


6.310 3 
2.810 2 
9.610 3 


2.410- 5 
6.810- 1 
1.210 3 


35 M , fast wind. Model from Garcfa-Segura et al. 1996b 
O-star 4.510 6 5.610~ 7 1.1 10 36 1.410" 3 
RSG 2.310 s 8. HO" 5 1.410 35 6.110 3 
W-R 1.910 s 2.410- 5 3.410 37 5.010" 2 


1.610 s0 
1.110 48 
2.010 s0 


6.310 3 
1.410 3 
9.610 3 


2.410- 5 
1.410 -1 
1.210~ 3 


60 M . Model from Langer et al. 1994 
O-star 1.7 1 6 4.7 1 0~ 6 1.010 37 
pre-LBV 1.510 6 1.5 1 5 2.710 37 
LBV 1.210 4 6.710" 4 2.610 37 
W-R 6.510 s 2.410~ s 3.010 37 


1.210- 2 
3.610~ 2 
2.410- 1 
4.810-2 


5.410 s0 
1.410 s1 
9.910 48 
6.210 s0 


2.1 10 4 
5.810 4 
2.810 3 
3.1 10 4 


2.5 10~ 4 
9.210" 4 
2.410" 1 
2.010- 3 
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Figure captions 

Fig. 1. CO and 13 CO integrated emission from 50 to 56 km s _1 toward NGC 2359. a) J-band picture of the optical nebula, 
overlaid with the CO J = 2 — » 1 emission from Rizzo et al. (2001a). The rectangle indicates the region mapped in this paper with 
higher angular resolution, b), c) and d) Molecular transitions indicated at the top left of each map, with the beams plotted in the 
bottom left corner. Angular offsets in Figs, b), c) and d) are referred to (RA, Dec.) 20 oo = (07 h 18 m 35 s :7, -13°16'37'.'0). First level 
and spacing in these figures are 4.4, 0.8 and 4 Kkm s _1 , respectively. 

Fig. 2. Comparison of the CO emission in two velocity channels, a) CO J = 1 — > brightness temperature at 55.2 km s _1 (grey 
contours), overlaid with the emission at 50. 1 km s - ' (black contours). Lower contour and step are 0.3 K (grey) and 0.16 K (black), 
b) The same as a) for the CO J = 2 — > 1 brightness temperature. Lower contour and step are 0. 15 K (grey) and 0.8 K (black). We 
note in both maps a clear shift in the emission, separated by ~ 20". 

Fig. 3. A position-velocity strip in the CO J = 2 — ♦ 1 emission, a) CO J = 2 — » 1 emission map integrated between 52 and 
55 km s _1 (on greyscale, from 2 to 22 Kkm s _1 ) and between 48 and 51 km s _1 (contours, 1 to 5 times 0.2 Kkm s _1 ). The 
straight line indicates the location of the strip. The filled squares along the line show the individual positions observed with 
higher sensitivity (see Sect. 3.2 and Fig. 4). The three kinematical components described in Sect. 4 are distinguished, since the 
greyscale represents the low velocity component, the contours sketches the intermediate velocity component and the high velocity 
component is detected just in position E. b) Position-velocity strip in the velocity range of interest. Angular distance is measured 
from southeast to northwest. The bottom small triangles indicate the approximate location of the individual positions marked by 
squares in a). Contours are 1 to 10 times 0.6 K. In general, we do not see a velocity variation greater than 1 km s _1 along the 
slice. However, a striking line broadening toward the inner part of the nebula appears, in close coincidence with the thin layer at 
48-50 km s _1 depicted in Fig. 2. 

Fig. 4. High-sensitivity, individual spectra of CO isotopes along the strip shown in Fig. 3. Every column corresponds to the 
spectral line written at the top, while every row corresponds to the selected positions. The assymetries in positions C, D and E 
indicates the presence of different velocity components. It is striking the very broad emission of the CO J = 2 — » 1 line in position 
E. 13 CO J = 2 — > 1 and C ls O J = 1 — » are detected for the first time in a W-R nebula. Temperature scale is T^ (K). Since 
different scales are used, auxiliary lines at 1 K (for both CO lines) and 0.2 K (for the 13 CO lines) were also plotted. 

Fig. 5. CO J = 2 — > 1 to J = 1 — > line ratio (R eX c) in the southern part of NGC 2359. The velocity ranges are indicated at the 
left of each figure, a) "low velocity component", b) "intermediate velocity component". In general, R exc does not show no large 
changes in both velocity intervals, except toward the inner nebula, in panel a), where R exc is at least a factor of six larger than the 
mean value. 

Fig. 6. CO to 13 CO J = 1 -» line ratio (R isot ) in the southern part of NGC 2359. The velocity range of integration is (52, 55) 
km s _1 . A systemic change (almost by one order of magnitude) from east to west is noted all over the map. This might indicate a 
different distribution of CO and 13 CO, variation of abundance or opacity effects in the area. 

Fig. 7. Variation of some physical parameters for the low {Ivc; filled squares) and intermediate (ivc; filled triangles) velocity 
components, as a function of the distance to HD 56925. a) Distribution of the density (in logarithmic units), b) 13 CO column 
density. Abcissae are projected distance to the W-R star by assuming a distance of 5 kpc. Open symbols represent upper limits. 



